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Abstract 

Widely studied quasiparticle models for quark gluon plasma is revisited here to understand the 
statistical mechanics and thermodynamics of the system. We investigate the statistical mechanics and 
thermodynamics inconsistencies involved in these models and their consequences in the observables. 
Quasiparticle model is a phenomenological model with few parameters and by adjusting them all models 
fit the results of lattice gauge simulation of gluon plasma 1 . However, after fixing 2 of the 3 parameters of 
the model by physical arguments, only one quasiparticle model, which is consistent with both statistical 
mechanics and thermodynamics, fits the Bielefeld lattice data pQ. The same model also fits the recent 
lattice results of Wuppertal-Budapest group [SJ, which deals with precision SU(3) thermodynamics for a 
large temperature range, reasonably well. 
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1 Introduction 



Quark gluon plasma (QGP) is a non-ideal system and hence the formulation of statistical mechanics (SM)and 
thermodynamics (TD) is a challenging one. SM and TD of ideal system is straight forward, but the equilib- 
rium studies of non-ideal system is still underdeveloped because of the difficulties due to mutual interaction 
among the constituents. Of course, there are theories like classical and quantum Mayer cluster expansion 
[3], but they are not fully successful for QGP because of the strong interaction [3]. Mayer's cluster expan- 
sion is valid for weakly coupled system, whereas QGP seems to be strongly coupled [S]. Therefore, one 
goes for phenomenological models like strongly coupled quark gluon plasma [SJ [3], quasi-particle models 
0IS![lOl[TlJ[T2[T3[T![T3[^ etc. with few fitting parameters. Model results are fitted 

with the results of the simulation of lattice gauge theory (LGT) for QGP [TJ [23 Here we study and 
compare different quasiparticle models and discuss SM and TD consistencies of each models. 

Next section we discuss the quasiparticle model which was hrst introduced. The modification of this 
model by eliminating the TD inconsistency is presented in Section III. Another alternative way without TD 
and SM inconsistencies is presented in Section IV. Discussions and conclusions is presented in Section V. 



Quasi-particle model for QGP was first introduced by Goloviznin and Satz [7], and then Peshier et. al. for 
gluon plasma [5] . Following the recent article by Castorina, Miller and Satz [2H H3] , one starts from the 
expression for pressure, 



where df = and T is the temperature. K2{x) is the modified Bcssel function of order 2. in is the thermal 
mass which is a function of temperature. Basic idea behind the quasiparticle model is that the effects of mu- 
tual interaction among the constituents in nonideal plasma is taken in the form of effective mass or thermal 
mass and treat the system as ideal. There are different phenomenological models for thermal masses and 
here we use the recent model, [22] , 



2 Quasiparticle Model - I (QPM-I): 





(i) 



l 



where t 




'I"ik: structure of this m(T) is discussed, in detail, in Ref. 2% which involves the main properties of phase 
transition near critical temperature and proper high temperature limit for large temperature. 

All other thermodynamic quantities like energy density, speed of sound, etc. are derived from pressure 
using thermodynamic relations. For e.g., the energy density, e, is 



Let us call this model QPM-I. 

However, above expression for energy density leads to statistical mechanics inconsistency. We know that 
when we treat an equilibrium system as canonical ensemble (CE) or grand canonical ensemble (GCE) the en- 
ergy of the system is not constant and fluctuates and we define average energy [3] and hence energy density as 



Therefore, we see that e ^ £sm and hence SM inconsistency. This model violates the standard SM of CE 
or GCE system. 

But, since m(T) has few parameters, we can tune them to fit LGT data for pressure and them obtain 
other TD quantities, using TD relations. Hence it is TD consistent. Thus we see that this model is TD 
consistent, but SM inconsistent. 



In this model, QPM-II, we follow Peshier et. al. [Sj where one starts from pressure, just like the earlier model, 



(3) 




(4) 



3 Quasiparticle Model - II (QPM-II): 




(5) 



However, the energy density e is 




(6) 



2 



5 
4.5 

4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 




— r - 



-i 1 - 



— H -x-- 



1.5 



2.5 3 
T/Tc 



3.5 



4.5 



Figure 1: Plot of e/T A and P/T 4 as a function of T/T c for gluon plasma from 3 quasiparticle models QPM-I, 
QPM-II, QPM-III and compared with LGT results of Ref. [I] (symbols). 

which is different from QPM-I. Hence it is SM consistent. But the drawback is that it is TD inconsistent as 
pointed out by various authors O EH E] ■ 

To solve this problem, borrowing idea from the bag model of hadron spectroscopy, one introduces bag 
pressure term B(T) in the expression for pressure as, 



P = d f T*f*A' Ka & ) 



l=i 



I 4 



B(T) 



(7) 



and energy density, 



oo 1 



l=l 



B(T) 



(8) 



This B(T) was interpreted as temperature dependent bag pressure, vacuum pressure, residual interaction 
energy, etc. Then one exploit B(T) to solve the TD inconsistency problem by demanding that pressure and 
energy density, Eqs. (|7J| and (JSJ), be TD consistent. This leads to TD consistent relation, 



dB 
dT 



l=l 



l=l 

and hence BIT) is related to m(T). Thus the QPM-II is now TD and SM consistent. 



(9) 
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Figure 2: Plot of e/T 4 andP/T 4 as a function of T/T c for gluon plasma from 3 quasiparticle models QPM-I, 
QPM-II, QPM-III at Boltzmann limit and compared with LGT results of Ref. [I] (symbols). 

4 Quasiparticle Model - III (QPM-III): 

Here we follow yet another approach, QPM-III, which follow straight forward from standard SM. Once the 
interaction effects of nonideal system is absorbed in thermal mass and treat the resultant system as ideal 
system in CE or GCE. We know that in CE or GCE, one defines average energy [3] and hence we start from 
energy density, 



1 = 1 



(10) 



and hence SM consistent. Pressure is obtained from the TD relation, 



dT 



(11) 



on integration and thus it is TD consistent also. Therefore, this model is both SM and TD consistent and 
no need of B(T). Of course, one can include B(T) in this model also, but then B(T) is free function of T. 
As we have shown in Ref. 24 this model gives QPM-II with the help of an extra constraint equation, Eq. 
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QPM-I 
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QPM-II 
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QPM-III 
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Table 1 



5 Discussions and conclusions: 

Thus we saw that QPM-I is SM inconsistent, but by exploiting the phenomenological parameters we can 
fit the LGT data pQ of gluon plasma for pressure. Values of parameters are tabulated in Table-I. But it is 
interesting that even after the violation of SM, by adjusting three parameters one can get reasonably good 
fit as shown in Fig. 1. Using the 3 fitting parameters other TD quantities like e is evaluated from the TD 
relations and also plotted in Fig. 1 and compared with LGT data and fits well as it should. 

Next we discuss the results of QPM-II. It is SM consistent and also TD consistent by TD consistent 
relation, and again reasonably good fit to LGT data for pressure and energy density are obtained by adjusting 
3 parameters which are tabulated in Table-I. Plots are given in Fig. 1. 

Results of QPM-III, which is both SM and TD consistent, is also plotted in Fig. 1, by fitting energy 
density with LGT data by adjusting 3 parameters. Values of fitted parameters are given in Table -I. Again 
reasonably good fits are obtained. 

It is interesting to note that even the single term in the summation for the expression of thermodynamic 
quantities, pressure or energy density, in all three models, also fits the LGT data by different sets of fitting 
parameters given in Table-I (right). These results are equivalent to the Boltzmann limit discussed in Ref. 
[22] . Plots are compared in Fig. 2. 

In the Ref. |22) the reason for the 3 parameters in m(T) is discussed, c is related to the critical exponent 
of phase transition and is estimated to be 0.41 for SU(2). Parameter b is introduced based on the fact 
that thermal mass is supposed to be proportional to T as T — > oo. We may even demand that at high 
temperature m(T) may goes as \J\^>p where ui p is the plasma frequency which follows from perturbativc 
results of quantumchromodynamics (QCD). Therefore, we may take 

m(T) = -^—+b(t)t, (12) 

where b(t)/T c = ^2Tra s (t), which is not constant as in Ref. [22], but running constant because of a s (T). We 
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Figure 3: Plot of e/T 4 and P/T 4 as a function of T/T c for gluon plasma from our model QPM-III with 
running b(t) and compared with LGT results of Ref. PQ (symbols). 



used the QCD perturbative result, uj 2 p = 47ra s (T)T 2 /3. We may take 2-loop perturbative coupling constant, 



()7T 



. 3(153 - 1971/) ln(2 ln(T/A T )) 



(13) 



4tt (33-2n/)ln(T/A T ) \ (33 -2n f ) 2 ln(T/A T ) 
where At = 0.14 T c from the lattice studies Ref. [T] with which we compare our results. Assuming that 
SU(3) also have the similar phase transition as SU(2), we may take c = 0.41 for gluon plasma. Now we have 
only one free parameter, a, and by adjusting it to fit to LGT data using our model QPM-III, we again get 
good fit as plotted in Fig. 3 for a = 0.45. Advantage here is that 2 of 3 parameters of m(T) is fixed by 
physical arguments. At present we don't know to estimate the value of a. We also tried to fit LGT data 
using QPM-I and QPM-II with the same At = 0.14 and c = 0.41, but by adjusting a. However, both of 
these models fail to get good fit to LGT data. 

It is further interesting to note that the same model, QPM-III, also fits the recent lattice results of 
Wuppertal-Budapest group [2 , which deals with precision SU(3) thermodynamics for a large temperature 
range, reasonably well as shown in Fig. 4. Here also a = 0.45 and c = 0.41, but ot s (T) is the 3- loop running 
coupling constant used in the same reference [2] . Our result is less than 2 % lower than that of lattice for 
large temperature as can be seen from Fig. 4. There is no much difference even if we use earlier a s (T), 
Eq.flTU), with the same A T = 0.14 T c . 
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Figure 4: Plot of e/T 4 and P/T 4 as a function of T/T c for gluon plasma from our model QPM-III with 
running b(t) and compared with LGT results of Ref. [2[ (symbols). 

In conclusion, all three quasiparticle models [22] EH EQ, considered here, fit the LGT data on gluon 
plasma [I] reasonably well by adjusting three parameters of a phenomcnological thermal mass proposed 
recently by Castorina, Miller and Satz [22] ■ However, only one model [IT], QPM-III, which is both statistical 
mechanics and thermodynamics consistent, fits the LGT data after fixing 2 of 3 parameters by physical 
arguments and adjusting the remaining one free parameter. It is remarkable that the same model, QPM-III, 
which fits Bielefeld lattice data |JJ , also gives good fit to the recent lattice results of Wuppertal-Budapest 
group [2] which includes data for larger range of temperature (> 5T C ). 
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